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ABSTRACT 

Context. A longstanding problem in astrochemistry is how molecules can be maintained 
in the gas phase in dense inter- and circumstellar regions at temperatures well below 
their thermal desorption values. Photodesorption is a non-thermal desorption mecha- 
nism, which may explain the small amounts of observed cold gas in cloud cores and 
disk mid-planes. 

Aims. This study aims to determine the UV photodesorption yields and to constrain 
the photodesorption mechanisms of three astrochemically relevant ices: CO, N2 and 
CO2. In addition, the possibility of co-desorption in mixed and layered CO:N2 ices is 
explored. 

Methods. The UV photodesorption of ices is studied experimentally under ultra high 
vacuum conditions and at astrochemically relevant temperatures (15-60 K) using a 
hydrogen discharge lamp (7-10.5 eV). The ice desorption is monitored by reflection 
absorption infrared spectroscopy of the ice and simultaneous mass spectrometry of the 
desorbed molecules. 

Results. Both the UV photodesorption yield per incident photon and the photodesorp- 
tion mechanism are highly molecule specific. The CO photodesorbs without dissociation 
from the surface layer of the ice, and N2, which lacks a dipole allowed electronic tran- 
sition in the wavelength range of the lamp, has a photodesorption yield that is more 
than an order of magnitude lower. This yield increases significantly due to co-desorption 
when N2 is mixed in with, or layered on top of, CO ice. CO2 photodesorbs through 
dissociation and subsequent recombination from the top 10 layers of the ice. At low tem- 
peratures (15 - 18 K), the derived photodesorption yields are 2.7(±1.3) x 10 -3 and < 
2 x 10 -4 molecules photon -1 for pure CO and N2, respectively. The CO2 photodesorp- 
tion yield is 1.2(±0.7) x 10 -3 x(l- e - (a:/2 - 9C±1 - 1)3 )+l.l(±0.7) x 10 -3 x (l-e" (a:/4 ' 6(±2 - 2) )) 
molecules photon - , where x is the ice thickness in monolayers and the two parts of 
the expression represent a CO2 and a CO photodesorption pathway, respectively. At 
higher temperatures, the CO ice photodesorption yield decreases, while that of CO2 
increases. 
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1. Introduction 



In dark clouds molecules and atoms collide with and stick to cold submicron-sized dust 



particles, resulting in icy mantles (jLeger et al 



1985 



Boogert fc Ehrenfreund 



20041 ). The 



ices are subsequently processed by atom or ligh t interactions to form more complex species 



( Tielens &: Hagen 



1982 



Watanabe et al 



2003 : ?). Observations show that H 2 0, CO and 



CO2 are the main ice constituents, with abundances up to 10~ 4 with respect to the total 
hydrogen de nsity. These molecules are key constituents in the formation of more com- 
plex species (jTielens fc Charnlevlll9971 ). and their partitioning between the grain and gas 
phase therefore stro ngly affects the chemical evolution in star- and planet-forming regions 



( van Dishoeck 



2006) 



Whether formed on the grains or frozen out from the gas phase, chemical models of cloud 
cores show that all molecules except for H2 are remove d from the gas phase wit hin ~ 10 9 /rtH 
years, where «h is the total hydrogen number density ( Willacy fc Millarlll998 ). For a typical 
cloud core density of 10 4 cm -3 , this time scale is much shorter than the estimated age of 
such regions and thus molecules like CO and CO2 should be completely frozen ou t. Yet 



gas-phase molecules, like CO, are detected in these clouds (jBergin et al 



CO gas is also d etected in the midplanes of protoplanetary disks (jDartois et al 



Pietu et al 



2001 



2002). Cold 



2003 



20071 ) where the densities are higher and the freeze-out time scales are even 



shorter, suggesting the existence of either efficient non -thermal desorption or an efficient 

(|2008f) have detected cold HCO^ , tracing 



Sakai et al 



mixing process in the disks. Similarly 
gas phase CO2, toward the embedded low-mass protostar IRAS 04368+2557 in L 1527 also 
referred to as L 1527 IRS. From the high column density and the thin line profile they 
conclude that the observed CO2 cannot originate from thermal evaporation of ices in the hot 
inner regions of the envelope. They instead suggest gas phase formation of CO2 to explain 
their observations, but do not consid er non-thermal deso rption in the cold envelope as an 

(jl999l ) toward several small translucent 



Turner et al 



alternative. HCOj is also detected by 
molecular clouds. They conclude that the observed HCOj can only form through gas phase 
chemistry for very specific C/O ratios and time spans and that the source of gas phase 
HCOj may instead be desorbed CO2 ice. Both the CO and CO2 observations may thus 
be explained by non-thermal desorption of ices, but this has not been quantified to date. 

In dense clouds and in outer disks and disk midplanes, desorption must occur non- 
thermally since the grain temperature is low enough, around 10 K, that thermal desorp- 
tion is negligible. Suggested non-thermal desorption pathways include photon and cosmic 



ray induced processes and desorption following the release of chemical energy (jShen et al 



2004; 



Roberts et al 



2007|). The importance of these processes depend both on the intrinsic 



desorption yields and on the local environment, especially the UV and cosmic ray fluxes. 
External UV photons from the interstellar radiation field can penetrate into the outer re- 
gions of dense clouds and disks and this UV field may be enhanced by orders of magnitude 
in disks through irradiation by the young star. In addition to direct interaction with ices, 
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cosmic rays and X-rays also produce a UV field inside of the clouds through interaction 
with H2. 

UV photodesorption is therefore possible in most dense astrophysical environments, 
but it has been proposed as an important desorption pathway of ices mainly in protoplan- 
etary dis ks and other astrophysical regions with dense clumps of material and excess UV 



photons l| Willacy fc L anger 



2000 



Dominik et al 



2005). There is however a lack of exper- 



imentally determined photodesorption yields for most astrophysically relevant molecules. 
This has prevented progress in the field and in most models UV photodesorption is simply 
neglected. Recently we show ed that CO photod esorption is an efficient process with a yield 
of 3(±1) x 10~ 3 photon -1 Joberg et all 120071). T his is of the same order as H 2 pho- 
todesorption, investigated by 



West lev et al 



(|1995al lbl). though the dependence of the H2O 
yield on different parameters remains unclear. The p hotodesorption of H2 O and benzene 
in a H2O dominated ice has also been investigated by [Thrower et al" ( 2008 ) who only find 
substrate and matrix mediated desorption processes. 

In this study we determine the photodesorption yield of CO2 and its dependence on ice 
thickness, temperature, morphology, UV flux and integrated UV flux or fluence as well as 
UV irradiation time. In addition, we extend the previously reported investigation of CO 
and N2 photodesorption to include different temperatures and ice morphologies. From the 
deduced yield dependencies we constrain the different desorption mechanisms and discuss 
the astrophysical implications. 



2. Experiments and their analysis 

2.1. Experimental details 



The experimental set-up (CRYOPAD) is described in detail by iFuchs et alj (|2006l ) and 



Oberg et al 



(2007). The set-up allows simultaneous detection of molecules in the gas phase 
by quadrupole mass spectrometry (QMS) and in the ice by reflection absorption infrared 
spectroscopy (RAIRS), with an angle of incidence of 84°, using a Fourier transform infrared 
(FTIR) spectrometer. The FTIR covers 1200 - 4000 cm -1 with a spectral resolution of 0.5- 
1 cm -1 . 

In the experiments presented here, thin ices of 2.1-16.5 monolayers (ML) are grown 
with monolayer precision under ultra-high vacuum conditions (P ~ 10~ 10 mbar) at 15 - 
60 K on a gold substrate that is mounted on a He cryostat. All experiments are conducted 
with the isotopologues 13 CO and 13 C 18 (Cambridge Isotope Laboratories 99% and 97% 
purity, respectively), 15 N 2 (Campro Scientific 98% purity), and 13 C02 (Indugas 99% pu- 
rity) and 13 C 18 02 (ICON Isotopes 96% purity) to avoid contributions from atmospheric 
contaminations as well as to be able to separate CO and N2 mass spectrometrically with the 
QMS. Test experiments show that the isotopologue choices do not affect the experimental 
outcomes for any of the ices. 

Within the uncertainties of the experiment, we also find that there is no difference in 
the photodesorption yield of 6.5 ML CO2 ice deposited on top of another 7 ML CO2 ice 
(of a different isotopologue), or 7 ML CO2 on top of 10 ML of H2O ice, compared with 
6-7 ML CO2 ice deposited directly onto the gold substrate. Since the character of the 
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substrate seems to have no influence on the photodesorption process, all other experiments 
are carried out with CO2 ices directly on top of the gold substrate. 

The ice films are irradiated at normal or 45° incidence with UV light from a broadband 
hydrogen microwave discharge lamp, which peaks a round Ly a at 121 nm and covers 115- 



170 nm or 7-10.5 eV (Munoz Caro fc Schuttc 



2003T ). All photodesorption experiments are 



performed in the same experimental chamber and the different UV angles of incidence 
are obtained by rotating the gold substrate. The lamp emission resembles the spectral 
distribution of the UV interstellar radiation field that impinges externally on all clouds. 
It is also consistent with the UV radiation produced locally inside clouds by the decay of 
electronic states of H2, following excitat ion by energetic e l ectro ns resulting from cosmic-ray 



induced ionization of hydrogen, see e.g. 



Sternberg et al 



(ll987D . 



The 45° and the normal incidence irradiation settings produce the same experimental 
results, except for a reduced photon flux on the ice in the 45° setting due to geometry. 
The lamp UV flux is varied between 1.1 and 8.3 xlO 13 photons cm -2 s _1 in the different 
experiments. The UV flux is monitored during all experiments using the photoelectric effect 
in a thin gold wire in front of the lamp. Before the start of the experimental run, this gold 
wire current was calibrated to an absolute flux in a separate set-up by simultaneously 
measuring the flux with a NIST calibrated silicon diode and the current induced in the 
gold wire. During the photodesorption experiments the flux onto the ice surface is also 
estimated by measuring the CO2 photodissociation cross-section several times during the 
experimental run, at both normal an d 45° incidence, and comparing our derived cross 

( 20031 ) . This was deemed necessary 



sections with the calibrated values in 



Cottin et al 



since the calibration measurements were carried out with normal incidence in the separate 
set-up, while most experiments used a 45° incidence angle. To prevent photodesorption 
during these measurements, the CO2 ice is covered with an inert ice layer. We find that in 
the normal incidence setting the resulting flux using this actiometry method deviates by a 
factor 0.9-1.4 from the photodiode-calibrated gold-wire results. 

Tables Q] and [2] summarize the experiments in this study. In the CO experiments, the 
temperat ure is varied betwe en 15 and 27 K, which is close to its thermal desorption tem- 
perature (jOberg et al.ll2005T ). This complements the previous CO photodesorption experi- 
ments, which investigated the d ependences of photodesorption on lamp flux and ice thick- 
ness at 15 K ( Oberg et al.ll2007 ). In three additional experiments the photodesorption yield 
(or upper limits) of N 2 is determined, as well as the changes in CO and N 2 ice photodesorp- 
tion in a CO:N2 ice mixture and in a N2/CO layered ice at 16 K. In the CO2 experiments 
the temperature is set to 16 - 60 K, the irradiation flux to 1.1 — 8.3 x 10 13 photons cm~ 2 
s _1 and the ice thickness to 2.1 - 16.5 ML. 



2.2. Data analysis 

The UV induced ice loss yield during each CO and CO2 experiment is determined by RAIRS 
of the ice during irradiation. The intensity of the RAIRS profile is linearly correlated with 
the ice layer thickness of CO and CO2 ice up to ~5 ML, but the RAIRS profile can be 
used up to 20 ML for analysis as long as the non-linear growth above 5 ML is taken into 
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Table 1. Summary of CO and N2 experiments 



5 



Experiment 


Composition 


Temperature (K) 


Thickness (ML) 


Lamp flux (10 13 photons cm 2 s 1 ) 


1 


13 CO 


15 


4 


4.7 


2 


13 CO 


22 


3.5 


4.7 


3 


13 CO 


27 


3.5 


4.7 


4 


i3 CO a 


16 


3 


4.7 


5 


15 N 2 


16 


4 


4.7 


6 


13 CO: 15 N 2 mixed 


16 


4:4 


4.7 


7 


15 N 2 / 13 CO layered 


16 


1/4 


4.7 



"annealed ice (i.e. deposited at 27 K and then cooled down to 16 K) 




5 10 15 20 

Deposited ice / ML 

Fig. 1. The integrated absorbance of the 13 C 18 02 stretching band as a function of deposited 
ice. The fitted exponential function is used to correct for non-linear growth of the integrated 
absorbance above 5 ML. Below 5 ML the RAIRS absorbance and the ice thickness are 
linearly correlated within the experimental uncertainties. 

account (Fig. [1]) . The absolute loss yield in number of molecules lost per incident photon 
is calculated from the RAIRS intensity loss as a function of UV fluence. 

For this calculation it is vital to have good estimates of the CO and CO2 band strengths. 
Due to the fact that all ice measurements are done using RAIRS, the ice thickness cannot 
be estimated from previously determined ice transmission band strengths. Instead the 
band strength of one ice monolayer is calculated from the observed difference in isothermal 
desorption from multilayer coverages (constant rate) and monolayer coverages (decreasing 
rate) as shown in Fig. [5] For CO this was done at 31 K and for CO2 at 76 K. The 
integrated absorbance of 1 ML is estimated to within 40% from the RAIR spectra at this 
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Table 2. Summary of CO2 experiments 



Experiment 


Composition 


Temperature (K) 


Thickness (ML) 


Lamp flux (10 13 photons cm 2 s 1 ) 


1 


13 c 18 o 2 


16 


5.5 


4.7 


2 


13 c 18 o 2 


18 


2.1 


2.3 


3 


13 c 18 o 2 


18 


5.5 


2.3 


4 


13 c 18 o 2 


18 


5.6 


2.3 


5 


13 c 18 o 2 


18 


9.0 


2.3 


6 


13 C 18 2 


18 


16.5 


2.3 


7 


13 c 18 o 2 


18 


3.9 


1.1 


8 


13 c 18 o 2 


18 


4.7 


3.5 


9 


13 c 18 o 2 


20 


6.2 


8.3 


10 


13 c 18 o 2 


18 


6.5 


8.3 


11 


13 C 18q 2 a 


18 


7.0 


2.3 


12 


13 C 18q 2 a 


16 


4 


4.7 


13 


13 c 18 o 2 


30 


6.2 


2.3 


14 


13 c 18 o 2 


40 


5.8 


2.3 


15 


13 c 18 o 2 


50 


6.7 


2.3 


16 


13 c 18 o 2 


60 


3.3 


2.3 


17 


13 c 18 o 2 


60 


7.4 


2.3 


18 


13 c 18 o 2 


60 


5.8 


1.1 


19 


13 c 18 o 2 


60 


6.2 


8.3 


20 


13 c 18 o 2 


60 


11.0 


2.3 


21 


13 co 2 


16 


3.8 


4.7 


22 


13 CO/ 13 C 18 2 b 


18 


10/5 


2.3 


23 


N 2 / 13 C 18 2 6 


18 


20/5.4 


2.3 


24 


13 c 18 o 2 /co 2 b 


18 


6.5/7 


2.3 


25 


13 C 18 2 /H 2 b 


18 


7/10 


2.3 



"annealed ice (i.e. deposited at 60 K and then cooled down to 16 or 18 K) 'layered ices 



turning point, which in its turn is used to calculate a band strength relevant for RA1RS. 
The calculated band strengths are 0.07 and 0.55 cm -1 ML -1 for CO and CO2 at their 
respective desorption temperatures. At 18 K the bands are somewhat weaker at 0.06 and 
0.42 cm^ 1 ML^ 1 , respectively. These values are highly set-up specific and they depend on 
such experimental parameters as mirror settings and should not be used for other purposes. 
This technique is based on the assumption that the ice is quite fla t at the desorptio n 



temperature, which was confirmed by the previous CO experiments (jOberg et al 



2007). 



For CO, the deduced ice thicknesses agree well (with in 20%) with the theore tical values for 
our chosen deposition pressure and deposition time (jAttard fc Barnesll2004l ). For CO2, the 
measured ice thickness is ^30% lower than predicted, probably due to the fact that some of 
the CO2 freezes out on the heating shield rather than depositing onto the substrate. Using 
this method we find that the relative band strengths of CO and C O2 ice compare well 



(within 20%) with previously measured transmission band strengths (jHudgins et al 



1993 
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Fig. 2. The integrated absorbance of the 13 C 18 (diamonds) and 13 C 18 02 (crosses) stretch- 
ing bands as a function of time during isothermal desorption of ~ 5 ML ices at 31 and 
76 K, respectively. The full lines are used to distinguish between the constant desorption 
rate in the multilayer regime and the decreasing desorption rate once the monolayer regime 
is reached. 



Gerakines et al 



19961 ). To convert the band strengths from cm 1 ML 1 to cm molecule 1 
a surface density of ~10 15 molecules cm -2 is assumed. 

In the case of CO, there is no measurable photodissociation in this wavelength range 



and the measured photon-induced loss yiel d is the photodesorption yield (jGerakines et al 



1996; 


Cottin et al. 


2003; 


Oberg et al. 


2007) 



20071 ). Simultaneous QMS measurements of the des- 
orbed CO gas phase molecules allow for the calibration of the QMS signal to an absolute 
photodesorption yield. The calibrated QMS signal for CO is used to determine the fraction 
of the CO2 ice that photodesorbs as CO. It is also used, together with the measured relative 
sensitivities of the QMS filament to CO and N2, to determine the N2 photodesorption yield. 
The N2 photodesorption yield cannot be determined by RAIRS since N2 has no permanent 
dipole moment and thus no strong IR feature. 

In contrast to CO, CO2 has only dissociative transitions in the wavelength region of the 
lamp; a UV photon ab sorption dissociates CO 2 into CO + O with a quantum yield of up to 
98% in the gas phase (Slangcr fc Black 19781 ). Hence, UV irradiation induces chemistry as 
well as desorption ( Gerakines et al. 19961 ). To determine the total photodesorption yield, 
the ice loss due to desorption must be separated from the conversion of CO2 into other 
ice products. This is done by analysis of the RAIR spectra using two different methods: 
kinetic modeling and mass balance. The first method uses the different kinetics of surface 
processes, like desorption, and bulk processes, like ice photolysis, to distinguish between 
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the two. Surface desorption from a multilayer ice is expected to be a zeroth order process 
and the photodesorption yield, which is determined from the derivative of the ice thickness 
with fluence, should therefore be constant with UV fluence. Photolysis into other species 
occurs throughout the ice at equal yield, since the ices here are thin enough that optical 
depth effects can be ignored, and is consequently expected to be a first order process. 
The contributions of desorption and photolysis to the observed ice loss is then determined 
by fitting the observed ice thickness versus fluence to the sum of a linear function and 
an exponential decay function, corresponding to a zeroth and a first order reaction. This 
method has the advantage that it is not dependent on identifying the photolysis products 
of CO2 and it is used to derive photodesorption yields whenever the zeroth and first order 
curves are separable. This is mainly the case for the high temperature and high fluence 
experiments. 

The mass balance analysis method compares the total ice loss with the simultaneous 
formation of other species in the ice; the final photodesorption yield is then defined as 
the loss yield of the original ice minus the formation yield of other carbon-bearing ice 
species. This is the only method that works for ices that are exposed to low fluences, where 
the contributions from the zeroth and first order processes cannot be s eparated. In the 



CO2 experiments CO, CO3 and O3 are the expected reaction products (jGerakines et al 



1996), with CO dominating. The photodesorption yield is then the CO2 loss yield minus 
the formation yield of CO and CO3, though as seen below the contribution from CO3 
is negligible. This method depends on the relative infrared band strengths of the formed 
molecules and is thus only accurate in the temperature range where the CO band strengths 
are measured, i.e. <30 K. As seen below, these two methods agree very well in the few 
cases where both methods are used to derive photodesorption yields. 

The simultaneous mass spectrometry of gas phase molecules during irradiation reveals 
the nature of the desorbed species. This is limited by the fact that less volatile molecules 
(e.g. CO2) adsorb onto the heating shield and other semi-cold surfaces inside the experiment 
before reaching the mass spectrometer. In the case of CO2, only the fractions of the ice 
that desorb as CO and other volatile species are detected by the QMS. At temperatures 
above 30 K the conversion factor between QMS detected and desorbed CO changes due to 
a decrease in cryopumping of CO, which is accounted for when deriving the CO-from-C02 
photodesorption yield. For both CO and CO2, re-condensation onto the actual ice sample 
after desorption will play a negligible role given the small surface area of the sample and 
the resulting underestimate of the actual photodesorption will be substantially lower than 
other sources of inaccuracy. 

To summarize, the main sources of uncertainty in these experiments are the photon flux 
and ice thickness calibrations of ^30% and 40%, respectively. In addition, from repeated 
experiments, the CO and CO2 experimental results are found to vary with approximately 
10% and 25%, respectively. The uncertainty is greater for CO2 than for CO because of 
the extra steps in deriving the CO2 photodesorption yield. The total uncertainty is ^60% 
for the CO2 photodesorption yield and ~50% for the CO photodesorption yield. The N2 
photodesorption yield uncertainty is somewhat higher due to the larger uncertainty in the 
QMS measurements. 
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Fig. 3. The CO ice photodesorption yield as a function of ice deposition temperature 
between 15 and 27 K. The ices are photodesorbed at the deposition temperature except 
for the point marked 27/16 K, where the ice is deposited at 27 K and then cooled down 
and irradiated at 16 K. The plotted uncertainties are the relative uncertainties between 
different experiments. The uncertainty in the absolute photodesorption yield is 50%. 



3. Results 



3.1. CO and N 2 

The results f rom p hotodesorption experiments of pure CO ice at 15 K are reported in 



Oberg et al 



( 20071 ). The yield presented there is updated here using new ice thickness 
and lamp flux calibrations; at 15 K the CO photodesorption yield is (2.7±1.3) xl0~ 3 
CO molecules per incident photon, averaged over the wavelength range of the lamp. The 
corresponding photodesorption quantum efficie ncy per absorbed UV photon in the surface 



Mason et al 



layer is estimated using t he lamp spect r um fr om lMunoz Caro fc Schuttd (|2003T ) and the UV 
spectrum of CO ice from 



(2006). The UV ice absorption spectrum shows that 
the CO lines are resolved. The measured absorption spectrum J" UV A^ c d\ is not calibrated 
to a UV cross section so in our calculation it is assumed that the total UV absorption cross 



sectio n is the same in the ice as measured in the gas phase J UV af^dX ( Eidelsberg et al 



1992t ?). The fraction of incident photons that are absorbed in the top monolayer, Xabs, is 
calculated from 

pUV UV-lamp AiccdA rUV eas^ 
Xabs = J „ TTV A .TTV-,. m „ ■ >< .TTV A . >< Ns (1) 



I 



UV r UV-lamp 
'A 



dA 



by cross-correlating the UV-lamp spectrum with the absorption spectrum of CO ice, di- 



vided by the total UV flux, J 



UV r UV-lamp 

'a 



dA, and then multiplying with the cross section 
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Fig. 4. CO spectra at 15 K (solid), 22 K (dotted), 27 K (dashed) and at 16 K and at 16 K 
of a sample deposited at 27 K (dash dotted), all acquired before irradiation. The figure also 
shows a spectrum of the annealed ice acquired after a UV fluence of ~ 7 x 10 17 photons 
cm~ 2 (thin dash dotted). 



conversion factor and the amount of molecules in one monolayer, N s . The resulting ab- 
sorption fraction is 5.5 ± 0.2 x 10~ 3 , where the uncertainty reflects the error in the gas 
phase cross section. Comparison with our measured photodesorption yield results in an 
efficiency of 0.3-0.8 per absorbed photon at 15 K, including both the absorbance and the 
photodesorption uncertainties. 

Figure [3] shows that the CO photodesorption yield decreases with ice temperature such 
that it is almost a factor of three lower at 27 K compared to 15 K. Within this temperature 
range the CO photodesorption yield is empirically described as linearly dependent on tem- 
perature: 2.7 x 10~ 3 — (T — 15) x 1.7 x 10~ 4 molecules photon -1 , where T is temperature 
in K. An additional experiment, where the ice is deposited at 27 K and then cooled down 
to 16 K before irradiation, results in a similar desorption yield as when the ice is also 
desorbed at 27 K. In quantum efficiency terms this corresponds to 0.1-0.3 photodesorption 
events per absorbed photon in the top ice layer, at 27 K as well as for the annealed ice 
at 16 K. This indicates that the structure of the ice, rather than the temperature, affects 
the photodesorption yield. This is further supported by a change in RAIRS profile at 22 
and 27 K compared to that at 15 K (F ig. HI). Changes in t he CO spectra with temperature 



have been previously reported by e.g. iFuchs et al.l (|2006l ). These spectral profiles do not 



change visibly after a UV fluence of 7 x 10 17 photons cm -2 when the ices are irradiated 
at their deposition temperature (not shown) . Figure 2] also shows the spectral profile of 
the annealed ice before and after irradiation, which has not changed significantly with cool 
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Fig. 5. The CO and N2 desorption yields and upper limits in 4 ML pure CO and pure N2 
ices, a 4:4 ML mixed ice, and a N2/CO 1/4 ML layered ice, all at 16 K, except for the pure 
CO ice at 15 K. In the mixed and layered ices the CO desorption yields are not constant 
and 'early' and 'late' marks the CO desorption yield in the beginning of the experiment 
and after a photon fluence of 8.5 x 10 17 cnr 2 . 



down, and it has at most slightly shifted toward the 15 K ice spectral profile following 
irradiation. 

In Oberg et al 



(|2007l ) the N2 photodesorption yield is constrained to a factor of 10 
less than the CO yield at 15 K. With increased sensitivity of the mass spectrometer, N2 
photodesorption is now detected at a yield of 1.8 x 10~ 4 N 2 molecules photon -1 or a factor 
of 15 lower than the CO photodesorption yield. This value has a factor of two uncertainty, 
mainly due to the uncertainty in the conversion of the mass spectrometer signal into an 
ice desorption yield. This measured yield is real, but because of continuous freeze-out of 
^0.1 ML H2O ice per hour the N2 ice contains an H2O impurity. A typical experiment 
lasts 5-6 hours resulting in a maximum 12% contamination level (a significant fraction of 
the adsorbed H2O molecules photodesorbs themselves during the irradiation experiments). 
This probably affects the measured photodesorption yield due to co-desorption of ices and 
thus the measured yield should be used as a strict upper limit of pure N2 photodesorption. 

In two additional experiments a CO:N2 mixture 4:4 ML and a N2/CO layered 1/4 ML 
ice are irradiated at 16 K. Figure [5] shows that in the ice mixture experiment the N2 
desorption yield more than doubles compared to pure N2 ice, while the CO desorption 
yield decreases by a factor of 2-4 compared to pure CO ice during the experiment. After a 
fluence of ~ 2 x 10 17 photons cm -2 the CO desorption yield is 50% of the photodesorption 
yield of pure CO. This is expected in a 1:1 mixture, since only 50% of the surface is covered 
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Wave numbers / cm" 1 

Fig. 6. 13 C 18 2 ice at 18 K before (bottom) and after (top) a UV fiuence of 5 x 10 17 
photons cm~ 2 . Some of the original CO2 ice is photolyzed into CO, CO3 (y\ at 1953 cm -1 
and 21/4 in Fermi resonance with v\ at 1810 cm -1 ). 

with CO. With increasing fiuence the CO desorption yield decreases such that it is only 25% 
of the yield of pure CO photodesorption after 8.5 x 10 17 photons cm~ 2 . This can only be 
understood if the N2 molecules desorb with a lower yield than the CO molecules, resulting 
in a decreasing CO surface concentration with UV fiuence. In the layered experiment the 
CO photodesorption yield is initially below the detection limit. After a UV fiuence of 
8.5 x 10 17 photons cm -2 the yield increases to 25% of the pure CO photodesorption yield. 
The N2 mass spectrometry signal does not reach equilibrium in the layered experiment, 
but the desorption yield seems to be at a similar level as in the mixed ice. 



3.2. C0 2 



3.2.1. Derivation of the total photodesorption yield 



To use the mass balance method to calculate the CO2 photodesorption yield, it is necessary 
to constrain which molecules are formed during irradiation. Figure [75] shows the spectra of 
an 18 K, 9 ML thick CO2 ice before and after a UV fiuence of 5 x 10 17 photons cm~ 2 . The 
only formed molecules are CO and CO3, though O3 formation cannot be excluded since the 
strong z/3 18 03 band around 1040 cm" 1 is outside of the range of the detector. This is in 

(|1996[ ). who found CO, CO3 and small am ounts of O3 after 



Gcrakincs et al 



agreement with 

irradia t ing CO2 wi t h a si milar fiuence. The line pos i tions are taken from 



Brewer fc Wang 



0972J), 



Moll et al 



( 19661) and iGerakines fc Moord (|200ll ). The weak band around 1605 



cm 1 cannot be unambiguously assigned. The lack of other features in the spectra from 
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Fig. 7. Top: the derived 13 C 18 2 ice thickness (diamonds) in a 18 K, 20/5.4 ML 
N2/ 13 C 18 02 layered ice as a function of fluence, plotted together with the formed 13 C 18 
ice (triangles) and the calculated total ice thickness (stars) . The latter is constant with flu- 
ence within the experimental uncertainties. The bottom panel shows the decreasing total 
ice thickness with UV fluence in an 18 K, 3.9 ML, uncovered experiment. In these plots 
the error bars indicate the relative uncertainty in the integrated absorbance (converted to 
a ML scale) of the RAIRS features within each experiment. This is also the case for similar 
plots throughout the paper 
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Fig. 8. RAIR spectra of the 13 C 18 2 stretching band at 2280 cm' 1 , the 13 C 18 stretching 
band at 2045 cm -1 and the 13 C 18 C>3 v\ band at 1950 cm -1 acquired before irradiation of a 
11 ML 13 C 18 2 ice at 18 K and then after every 8.3 x 10 16 photons cm~ 2 . The absorbance 
of the CO2 band decreases with UV fluence due to photodesorption and photodissociation, 
while the CO band absorbance increases. Note the nearly constant 13 C 18 03 integrated 
absorbance after a fluence of ~8xl0 16 photons cm~ 2 . 



e.g. carbon-suboxides put strict upper limits on the formation of such molecules to a 
fraction of a percentage of the original CO2 ice. This is consistent with the experiments of 



Gerakines fc Moord (|200lh where no carbon-suboxides was detected after UV irradiation 



of pure CO2 ice. In addition Temperature Programmed Desorption (TPD) experiments 
following irradiation show that ~10% of the original ice is photolyzed into CO, ~1% into 
O2 or O3 and ~0.1% into C2. From this we infer that more than 99% of the carbon budget 
is bound up in CO2, CO and CO3 during the experiment. 

The CO2 and CO abundances during each experiment are calculated from their derived 
band strengths. The CO3 band strengths have not been measured, however, and can only be 
crudely estimated. This will introduce a large uncertainty into the mass balance calculations 
if CO3 is formed at significant abundances. Figure [7J shows the calculated CO2 and CO 
ice thicknesses as a function of fluence for a layered 20/5.4 ML N2/CO2 ice. The ice cover 
hinders desorption and the result is that 10% of the original CO2 is photolyzed into CO. 
The fact that the lost CO2 is perfectly compensated for by the formed CO ice shows that 
CO3 is not a significant photolysis product. The lack of photodesorption in the layered 
experiment (top Fig. [7]) is contrasted with the observed photodesorption of a 3.9 ML bare 
C0 2 ice (bottom Fig.[7J. 
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Fig. 9. The calculated layer thicknesses of an originally 6.5 ML CO2 ice at 18 K as a 
function of UV fluence (diamonds) plotted together with the formed CO layer thickness 
(triangles) and the calculated total ice thickness (stars). The CO2 ice loss is fitted as a sum 
of photolysis (an exponential function) and desorption (a linear function). The exponential 
part mirrors the CO formation within the fit uncertainties. 



The CO3 abundance is also estimated independently by employing its only likely for- 
mation path and the fact that CO3 reaches its final level within ~ 5 x 10 16 photons cm~ 2 
in all experiments (exemplified in Fig. [5] where the level does not change between 0.8 and 
5xl0 17 photons cm -2 ). CO3 is expected to form from CO2 + O, where the O originates 
from photolysis of another CO2 molecule into CO + O. With photodesorption hindered, 
the amount of CO3 then never exceeds the amount of CO in the ice. In the 20/5.4 N2/CO2 
experiment, the CO3 abundance reaches steady state when less than one percent of the 
CO2 is converted into CO, which puts a 1% upper limit on the formed CO3. This is small 
compared to the ice loss during CO2 photodesorption experiments, where typically 20% of 
the ice is lost. 

From these results the mass balance photodesorption yield is defined as the CO2 ice 
loss yield minus the CO formation yield. Figure [5] shows the CO2, CO and CO2+CO ice 
thicknesses in an originally 6.5 ML thick CO2 ice at 18 K as a function of UV fluence. 
Practically the photodesorption yield is derived from the slope of the total ice thickness as 
a function of fluence. This mass balance method of determining the photodesorption yield 
agrees very well with the yield determined through simultaneous kinetic modeling of bulk 
photolysis and photodesorption. In Figure [5] the CO2 ice thickness is fitted to a function of 
the form A(0) x e ~ A W/* + A(2) + A(3) x $ using the IDL script MPFIT, where $ is the 
fluence. The photodesorption yield is determined from A(3). The derived photodesorption 
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Fig. 10. Mass spectra acquired during irradiation of a 6.2 ML thick 13 C 18 02 ice at 20 and 
60 K with a flux of 8.3xl0 13 photons cm -2 s _1 . The spectra at 60 K have been divided by 
3.3 to account for the lower cryopumping of volatiles like CO and O2 at 60 K compared to 
at 20 K. In each case the ice is irradiated for 3 hours before acquisition to ensure that the 
photodesorption rate is stable. Each acquisition lasts 3 hours and consists of ~100 averaged 
spectra. In addition to photodesorbed ices there are some background CO (m/z=12, 16 
and 28), CO2 (m/z=44) and possibly some background H2O as well (m/z=18). 



yield is the same, within the uncertainties, to the yield derived from fitting a linear function 
to the CO2+CO ice thickness. This confirms the validity of both methods. 

Whichever method is used, the result is a linear coefficient, which gives a pho- 
todesorption yield in loss of ice monolayers per 10 17 UV photons for a 6.5 ML ice at 
18 K: 0.27 (ML)/2.1 (10 17 photons cm- 2 ) = 0.13 ML / (10 17 UV photons cm- 2 ). This 
is further converted into a photodesorption yield in CO2 molecules per UV photon (7- 
10.5 eV): Y pd = 0.13 x 10" 17 ML photon" ^m 2 x (lO 15 molecules cm" 2 /l ML) = 1.3 x 10" 3 
molecules photon" 1 . 



3.2.2. Desorption products 

The total photodesorption yield is well constrained from the RAIR spectroscopy of the ice 
during irradiation. The question remains in which form CO2 ice photodesorbs. Figure [TTJ] 
shows mass spectra acquired during UV irradiation of a 6.2 ML ice at the two temperature 
extremes, 20 and 60 K. The only visible desorption products are CO and O2, which puts 
strict upper limits on other potential volatile desorption products. Less volatile species 
like CO2 (m/z=49 for 13 C 18 02) cannot be excluded, however, since their cryopumping 
efficiencies are up to two orders of magnitude higher than for CO. As described in Section 
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Fig. 11. The detected CO photodesorption during irradiation of a 4 ML 13 C 18 ice at 16 K 
and 13 C 18 2 ices at 18 K (5.5 ML) and at 60 K (7.4 ML), as a function of UV fluence. 
The 60 K signal has been divided by 3.3 to account for the lower cryo-pumping at 60 K 
compared to 18 K. 



2.2 the measured CO QMS signal can be converted into a number of CO molecules desorbed 
per photon. This number is compared with the total CO2 photodesorption yield to quantify 
the amount of the CO2 ice that desorbs as CO. From QMS measurements during irradiation 
experiments, ~ 20 — 50% desorbs as the fragment CO (Fig. QT]) and at most 5% as 02- It 
is thus inferred that more than 50% of the desorbed ice comes off as less volatile species, 
most likely CO2. This is supported by the fact that the amount of formed CO3 ice is the 
same whether or not the ice is covered and therefore whether or not photodesorption is 
allowed. This makes it unlikely that CO3 is photodesorbing in the uncovered experiments. 
In addition no other C-bearing molecules are formed at significant abundances, which only 
leaves CO2 as a possible desorption product. 

Below we separate the total CO2 ice photodesorption yield (as measured with RAIRS) 
from the desorption of CO-from-C02 (measured with the QMS). The desorption of CO2 
molecules is taken to be the difference between the total CO2 ice desorption and the 
CO-from-C02 desorption. The total photodesorption quantum efficiency per absorbed 
UV photon in the surface l ayer is estimated to 0.4-1 using the lamp spectra from 



Muii oz Caro fc Schuttd (|2003l ) and the calibrated UV spectra of CO2 ice from 



Mason et al 



( 2006T ). As seen from the thickness dependence below, this efficiency decreases with depth 
into the ice. 



18 



Karin I. Oberg et al.: Photodesorption of ices I: CO, N2 and CO2 




5 10 15 

Ice thickness / ML 



5 10 15 

Ice thickness / ML 



0.005 
0.004 
0.003 
0.002 
0.001 
0.000 



~ 1 — 1 — 1 — 1 — 1 — 1 — r 



n — 1 — i - 

_ 11 ML 



i(— 6 ML 



-3 ML 



0.005 
0.004 
0.003 
0.002 
0.001 



~i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 

d) 



t t * 



20 40 60 

Temperature /K 



0.000 
80 



J 1 1 1 I 1 1 1 I 1 1 1 L 



2 4 6 8 

UV flux / 10 13 photon cm" 2 s" 1 



10 



Fig. 12. Total CO2 (crosses) and CO-from-C02 (squares) photodesorption yield depen- 
dences on different parameters. In a) CO2 ices of different thicknesses are irradiated with 
the same UV fluence of ~6xl0 17 photons cm~ 2 at 18 K and in b) at 60 K. Both the 
total CO2 and the CO-from-C02 desorption yields are fitted with functions of the form 
c x (1 — e^WO) a t low temperatures (solid lines), where x is the ice thickness and / an 
ice diffusion parameter. The two measurements of 5.5 and 5.6 ML ices in a) are from the 
beginning and the end of a two-month long experimental series. Panel c) shows that the 
total photodesorption yield is constant with temperature within a low temperature region 
(<40 K) and within the warmer region 40-60 K for ^3 ML (diamonds), 6 ML (crosses) and 
11 ML (triangles) thick ices, irradiated with fluences of ~6xl0 17 photons cm~ 2 . Finally 
panel d) demonstrates the independence of the total and CO-from-C02 photodesorption 
yields on the photon flux for 4-6.5 ML ices at 18 K. 



3.2.3. Yield dependences on experimental parameters 

Ice thickness. The total CO2 photodesorption yield at 18 K is thickness dependent up to 
several monolayer s (Fig. | 12^l. which is in contrast with the constant photodesorption yield 



Oberg et al 



(|2007T ). The total CO2 photodesorption yield increases from 



of CO reported in 

8 x 10~ 4 to 2.3 x 10 -3 molecules photon -1 when the ice is grown from 2 to 16 ML. 

Figure [T2"k also shows simple models fitted to the measured yields of both the CO- 
from-C02 and the total CO2 photodesorption. The desorption yield of the CO-from-C02 
is well described by 1.1(±0.2) x 10~ 3 X (1 — e - ^/ 4 - 6 ^ 2 - 2 ^), where x is the ice thickness in 
monolayers. The total CO2 desorption dependence on thickness is a sum of the desorbed 
C0 2 and CO molecules and is thus modeled as 1.2(±0.1) x 10~ 3 x (1 - e -(*/2.9(±i.i))) + 
1.1(±0.2) x 10~ 3 x (i-e^/ 4 - 6 ^ 2 - 2 ))). Both models are fitted using the IDL routine MPFIT. 
Here the uncertainties only reflect the calculated model errors. As discussed in Section 2.2, 
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Fig. 13. The top panel shows the measured CO2 layer thickness (diamonds) of an originally 
7.4 ML CO2 ice at 60 K as a function of UV fluence, plotted together with the formed CO 
layer thickness (triangles) and the calculated total ice thickness (stars). The CO2 ice loss is 
fitted as a sum of photolysis (an exponential function) and desorption (a linear function). 
The bottom panel shows an 18 K experiment that has been thermally annealed at 60 K 
prior to irradiation at 18 K. Note the lack of evidence of photodesorption from this ice. 
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the total uncertainty is 60%. From these expressions it is clear that more than 90% of the 
photodesorption events originate in the top 10 ML and 50% in the top 3 ML of the ice. 
They also expose the thickness dependence of the fraction of the CO2 desorbing as CO; 
between 2 and 16 ML this fraction grows from 20 to 45%. The origin of the model and 
the full significance of the different exponential constants for CO and CO2 desorption is 
further discussed below. 

The build-up of molecules in the 18 K CO2 ice is linearly dependent on the ice thickness 
as expected for a photodesorption yield that is low in comparison to the total ice thickness. 
In all the ices, ~10% of the original ice is converted to CO and frozen into the CO2 ice 
after a typical UV fluence of 5 x 10 17 photons cm -2 . For comparison ^10% of the original 
ice is desorbed after the same fluence in a typical 5 ML experiment. 

Figure [T2b shows that the CO2 photodesorption has a somewhat different thickness 
dependence at 60 K. The CO-from-C02 desorption dependence on thickness is indistin- 
guishable from a linear function and is fitted linearly. The total CO2 desorption is fitted 
well, but not uniquely, by 2.2(±0.2) x 10 -3 x (1 - e -(*/5.8(±i.2))) + 2.2(±0.3) x 10 -4 x x 
molecules photon -1 , where x is the ice thickness in monolayers. This formula indicates 
that photodesorption takes place deeper into the ice at higher temperatures and that the 
mean-free-path of CO becomes infinite. 

Temperature. At temperatures higher than 30 K the CO2 photodesorption is initially 
fluence dependent, which is not the case for colder ices. This is further discussed below; here 
the constant yield reached after a fluence of 2.0 x 10 17 photons cm -2 is used for comparison 
between the photodesorption yields at different temperatures. Figure H"2"b shows the total 
CO2 photodesorption yield dependence on temperature for different ice thicknesses. At 18 
and 30 K both the total and the CO-from-C02 photodesorption yields are indistinguishable 
within the experimental uncertainties. Between 30 and 40 K the photodesorption yield 
jumps and above 40 K the photodesorption yield is again independent of temperature. 

The build-up of CO molecules in the CO2 ice is also temperature dependent. Above 
30 K, the CO build-up is less than 50% of the build-up at lower temperatures. 

Photon flux. Figure[T2"H shows 4-6.5 ML CO2 ices irradiated with different photon fluxes 
at 18 K. Between 1.1 and 8.3 x 10 13 photons cm -2 s -1 the CO2 photodesorption yield in 
molecules photon -1 is constant within the experimental uncertainties. The CO desorbing 
from the CO2 ice is independent of the lamp flux as well. At 60 K the ice is irradiated at 
two different fluxes and, similarly to the colder ice, the photodesorption yield is constant 
(not shown). 

At 18 K, the produced CO ice increases throughout the experiment up to a fluence of ~ 
10 x 10 17 photons cm -2 . After 10 x 10 17 photons cm -2 the amount of CO ice reaches steady- 
state, which is only clearly visible in the experiment with the highest fluence. The observed 
independence of the CO-from-C02 photodesorption yield on CO ice content indicates that 
direct CO photodesorption from the formed CO ice plays a minor role during irradiation 
of the CO2 ice. 

Time and photon fluence. Figure [TBI shows the photodesorption of a 60 K, 7.4 ML CO2 
ice irradiated with 2.3 x 10 13 photons s -1 cm -2 . At 60 K there is no measurable photodes- 
orption during the first two hours or a fluence of 2 x 10 17 cm -2 . This can be compared 
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with the 4.0 ML CO2 ice at 18 K in Fig. [7] (bottom), where photodesorption starts within 
a fluence of 10 17 photons cm -2 . In the 60 K experiment, the total CO2 photodesorption 
yield jumps to 3.0 x 10 -3 molecules photon -1 after a fluence of 2 x 10 17 photons cm -2 
and remains constant for the remainder of the experiment. This photodesorption delay is 
observed in all 5-7 ML ices at 40-60 K and also at all thicknesses between 3-11 ML at 
60 K. The QMS measurements also show clear differences between the 18 K and the 60 K 
experiments (Fig. ITT]) . At 60 K, the CO QMS signal is lower than at 18 K during the first 
2xl0 17 photons cm -2 , corresponding to a yield ~lx 10 -4 molecules photon -1 . After the 
first 10 17 photons cm -2 the CO signal increases rapidly with fluence. 

To test whether this delay in photodesorption onset at 60 K is time or fluence dependent, 
a 5.8 ML ice is also irradiated at a 50% lower flux. In this experiment the photodesorption 
according to the RAIRS begins after the same photon fluence, which occurs after twice the 
amount of time compared with the experiment at a higher flux. 

Thermal annealing. In one experiment, a 7.0 ML thick ice is deposited at 60 K and 
subsequently cooled down to 18 K before starting the irradiation. The CO-from-C02 pho- 
todesorption yield is ~ 1 x 10 -4 photon -1 , while only an upper limit of 5 x 10 -4 photon -1 
is derived for the total CO2 photodesorption from the RAIRS measurements. This is signif- 
icantly lower compared to both unannealed experiments at 18 K and to 60 K experiments 
(Fig. [T3"]) . The CO-from-C02 QMS signal is similar to that from warm ices during the 
first 10 17 photons cm -2 . The CO ice build-up, as measured from RAIRS, is similar to the 
un-annealed ice at 18 K. 



4. Discussion 

4.1. CO and A/2 yields and mechanisms 
In 



Oberg et al 



(120071 ) we concluded that CO photodesorbs from the top one or two ice 
layers at 15 K. Recent theoretical work shows that CO only desorbs from the absolute 
surface layer (Takahashi & van Hemert, in prep.), which is supported by the new findings 
in this study. 

The experiments show that the photodesorption of CO depends on thermal annealing, 
such that annealing at a higher temperature results in a lower desorption yield. The an- 
nealing most likely results in a more compact ice with a smaller effective surface area as 
well as stronger bound molecules. These two effects then add up to decrease the quan- 
tum efficiency of the photodesorption process. The reason for the linear dependence with 
temperature is unclear and may be coincidental. This behavior cannot be extrapolated to 
lower temperatures, since a 15 K ice should be amorphous. 

The new photodesorption experi ments with N2 fur ther constrain the CO photodes- 
orption mechanism. Consistent with IQberg et al.l (|2OO70 there is no evidence of direct N2 



photodesorption. The increase in photodesorption yield of N2 when mixed with CO or 
grown in a monolayer on top indicates that ~5% of the UV photon absorptions of CO 
molecules result in the desorption of a neighboring molecule rather than the desorption 
of the originally excited molecule. The decreasing CO photodesorption in the mixed ice 
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and the lack of initial CO desorption in the layered experiment also confirm that CO only 
desorbs from the top ice layer. 



4.2. CO2 yield and mechanism 



In contrast to CO and H2O photodesorption (jAndersson et al.ll2006t ?, Takahaski & van 
Hemert, in prep.), the photodesorption mechanism of CO2 has not been theoretically ad- 
dressed yet. From the dependencies reported here, the mechanism may be constrained, 
however, and it is similar to that of H2O. CO2 photodissociates to CO+O, where the prod- 
ucts have excess energy. This is followed by both reactions in the ice to form the observed 
CO3, and recombination to C02- Some of these reaction products subsequently desorb. 



Westlev et al. 


1995a 


Oberg et al. 


2007 



2007, Oberg et al. accepted by ApJ). 



This independence is expected for single photon processes, but not for multi photon pro- 
cesses or desorption induced by excess heat from the lamp. This is consistent with the sug- 
gested mechanism of dissociation fragments and recombined molecules traveling through 
the ices before desorbing. 

The CO2 photodesorption yield is clearly thickness dependent, which is in contrast to 
the CO photodesorption from surfaces only. This difference can be understood from the 
fact that before a desorption event occurs, the CO2 molecule is dissociated into energetic 
products (whether concerned with the fragments or recombined molecules), which may 
travel through several monolayers of ice before they are quenched by the surrounding 
matrix. Assuming a homogeneous ice, the probability of a molecule with excess energy, from 
dissociation or recombination, reaching the ice surface and desorbing is only dependent on 
the excess energy, the diffusion properties of the molecule and the ice depth at which it 
receives its excess energy. The diffusion properties of a molecule are simplest described by 
the average distance the molecule travels through the ice before being stopped. Defining I 
as the average distance traveled by a molecule before quenched by the surrounding ice, the 
fraction of particles with excess energy that will move through a slab of ice of thickness 
x is expected to be proportional to e~ x ' 1 , assuming uniform photochemistry throughout 
the ice and that the direction the molecule travels is independent of ice depth. Integrating 
over the ice depth from to x, the total desorption of particles between and x is then 
proportional to 1 — e~ x l l . 

We find that this type of expression describes the photodesorption at low temperatures 
well, with an average travel distance or mean-free-path of 2.9 ML for the CO2 photodes- 
orption and 4.6 ML for the CO fragment desorption. The two different values for the 
mean-free-paths have large uncertainties and it is not clear that the difference is signifi- 
cant. The values are however consistent with the different sizes of CO and CO2, since the 
larger molecule is expected to be less mobile in the ice. The different mean-free-paths may 
also be partly due to different mechanisms, i.e. dissociation versus recombination, through 
which the CO fragment and the recombined CO2 molecule receive excess energy. 

The different mean-free-paths for CO and CO2 desorption from CO2 ice can also be 
consistent with a different desorption mechanism of CO2 molecules - momentum transfer 
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to a surface CO2 molecule from a smaller fragment originating in an underlying layer. This 
is observed in simulations by ? as an equally important photodesorption pathway for H2O 
compared to desorption of the recombined molecule. The ice thickness dependence of this 
kind of process remains to be explored, since it requires more complex models than the 
simple mean-free-path model presented here. 

Other desorption pathways of the CO fragment can however be ruled out from the ex- 
periments presented here. The desorbed CO molecules do not originate in photodesorption 
of CO ice, since the yield does not increase with an increasing fraction of CO in the ice - 
this fraction never reaches equilibrium during most low temperature experiments - while 
the CO QMS signal does. In addition, the underlying substrate seems to have no influence 
on the desorption yield; hence substrate mediated processes are excluded as well. 

At temperatures above the pure CO desorption temperature of ~30 K, the increased 
photodesorption yield is most likely due to the increased mobility of CO2 and CO in the 
ice. This is seen from the longer mean-free-path for desorbing CO2 molecules and the 
infinite mean- free-path of CO. The latter points to a very high mobility of CO in the ice 
at these temperatures, which results in that most of the produced CO thermally desorbs 
once formed through CO2 dissociation. It is also shown by the onset of O2 photodesorption 
(Fig. [TO]). The desorption of O2 is less than CO at any temperature (Section 3.2.1), but its 
mere presence shows that at 60 K oxygen atoms are very mobile and produced abundantly. 

A curious and not yet fully understood feature is the delayed onset of both CO and 
CO2 photodesorption from CO2 ices at 40-60 K. Since it is a fluence rather than a time 
effect it is probably due to a re-structuring of the ice induced by UV irradiation. This 
restructuring is indicated by a change in the infrared spectral profile (not shown), which 
seems more pronounced for the warmer ices compared to the colder ones. This proposed 
re-structuring may also be responsible for the apparent initial increase in CO2 in Fig. 1131 
Ice band strengths depend on ice structure (e.g. ?) and this increase in CO2 signal is most 
likely due to structural changes upon UV irradiation, modifying the strength of the CO2 
ice band, rather than more CO2 molecules adsorbing. Similar changes in band strength 
have been previously observed in CO and in H2O ices bombarded by ions (??) 

The importance of ice structure for the photodesorption yield is seen for CO as well, 
as discussed above. For CO2 this importance is further shown by the experiment on a 
thermally annealed ice, where the photodesorption yield is < 40% of the yield of an amor- 
phous ice. This is in contrast with the increased photodesorption yield for warm (40-60 K) 
ices, which are expected to have similar ice structures as the annealed ice. This new ice 
structure probably has a smaller effective ice surface compared to the 18 K amorphous ice, 
explaining the low photodesorption yield of the annealed ice. At temperatures above 40 K 
this decrease in surface does not affect the yield due to the high mobility of molecules in 
the ice. In summary, the CO2 photodesorption yield is both temperature and structure 
dependent. 
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4.3. Astrophysical significance 



4.3.1. CO and N 2 



Oberg et al 



(2007). 



The relevance of CO photodesorption at 15 K is discussed in detail in . 
The results here show that the previously reported yield is dependent on temperature 
and annealing between 15 and 27 K. The astrophysical significance of this dependence is 
probably limited, however. Because the desorption temperature and annealing temperature 
are similar, only a minor fraction of the pure CO ice will ever be annealed in astrophysical 
environments. Therefore the photodesorption yield derived at 15 K may be used for most 
purposes, also where a real temperature gradient exists. 

While a temperature gradient may not significantly reduce the CO photodesorption 
yield, a layer of N2 ice on top of the CO ice does. This layer may form either through later 
freeze-out or through selective photodesorption as shown in the experiments. Because of the 
peculiar photodesorption mechanism of CO ice, this study shows that a single monolayer 
of N2 ice initially reduces the CO photodesorption yield with more than 80%. After a 
UV fluence of 8.5 x 10 17 photons cm~ 2 (corresponding to ^300 years of irradiation at 
cloud edges and ~3 million years in cloud cores), the yield has increased to 25% of the 
pure CO ice photodesorption yield, indicating that 25% of the N2 layer is desorbed. This 
probably happens through CO co-desorption, since the photodesorption yield of pure N 2 
is too low to account for the observed desorption. N 2 co-desorption is also observed in a 
mixed CO:N 2 ice where N 2 photodesorbs at a yield of 3 x 10~ 4 photon -1 . This is still 
a very low photodesorption yield compared to e.g. the CO photodesorption yield and its 
astrophysical significance is doubtful. Other non-thermal desorption mechanisms such as 
cosmic ray induced spot heating will likely be more important. 



4.3.2. C0 2 

At low temperatures and for thick ices the C0 2 photodesorption yield is 2.3(±1.4) x 10~ 3 
photon -1 , which is almost identical to the CO photodesorption yield. Of the desorbed 
ice more than 50% desorbs in the form of C0 2 , while the remaining 20-50% desorbs as 
CO. At higher temperatures some C0 2 , up to 5%, also desorbs as 2 . This means that 
similar C0 2 and CO abundances are expected in regions where photodesorption dominates, 
assuming that the photodesorption yields of the two molecules are not much different for 
astrophysical ice morphologies compared to the pure ices studied here. This remains to be 
investigated, especially for the astroph ysically relevant cases of CO ice on top of C0 2 ice, 



and CO:C0 2 and C0 2 :H 2 mixtures (jPontoppidan et al 



2008). Based on our results we 



expect the recombined and thus energetic C0 2 molecule to penetrate at least as many CO 
ice layers during photodesorption as the 10 ML of C0 2 ice observed in this study, since CO 
forms a more loosely bound ice compared to C0 2 . The behavior of C0 2 in a H 2 matrix 
is more difficult to predict. Hence, until such a laboratory study exists we recommend to 
use the equation presented here for C0 2 photodesorption, with the possible modification 
of taking into account that only a fraction of the ice is C0 2 ice. 
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In the regions in clouds and disks where ices begin to form, the ice thickness dependence 
of the CO2 photodesorption needs to be taken into the account. As soon as the grain 
is covered with less than 10 ML of CO2 ice, the CO2 photodesorption yield decreases 
considerably with ice thickness according to 1.2 x 1CT 3 x (1 — e - x / 2 - 9 ) + 1.1 x 10~ 3 x (1 — 
e~ x / 4 - 6 ) at 18 K, where x is the ice thickness and the two parts in the yield expression are 
due to CO2 and CO desorption, respectively, during the irradiation of a CO2 ice. 

Heated or thermally annealed, i.e. heated and sub sequently cooled down, C O2 ice is 
observed toward many high- and low- mass protostars f 3Pontoppidan et al.ll2008l ). Because 
of the irreversibility of the infrared spectroscopic signature of heated CO2 ice, the two 
cannot be easily distinguished. At temperatures higher than the pure CO sublimation line 
at 25-30 K, the C0 2 photodesorption yield increases to 2.2 x 10~ 3 X (1 - e^/ 5 ' 8 ) + 0.22 x 
10~ 3 x x molecules photon^ 1 for the CO2 and the CO-fragment desorption, respectively. 
This results in a CO2 yield increase of at most a factor of two. In addition, possible 
annealing may decrease the photodesorption yield to a value that is less than 40% of the 
cold amorphous CO2 ice yield. The uncertainty in the photodesorption yield is thus an 
additional factor of 2-3 in thermally processed regions. 



4.3.3. CO2 astrophysical model 



CO2 is not detected directly in the radio regime due to its lack of a permanent dipole 
moment. It is instead traced by HCO^~ toward e.g. the protostar L 1527 IRS. As a test 
case we use a simple model to investigate whether the recen tly observed HCOjj " toward L 



Sakai et al 



(2008) estimated 



1527 IRS may be explained by photodesorption of CO2 ice. 
the gas phase CO2 abundance to be > 2.9 x 10 -7 with respect to H2, by using a simple 
reaction scheme for the CO2 to HCOj chemistry a nd assuming t hat th e CO2 is extended 



Furlan et al 



(|2008h report a C0 2 ice 



over the beam size of the IRAM 30 m telescope, 
abundance of 5 x 10~ 6 with respect to H2 in L 1527 IRS, which is used in our model. This 
abundance is almost an order of magnitude lower than what is observed toward a large 
sample of low-mass protostars (jPontoppidan et al.ll2008l ) and may be underestimated. The 
effect of increasing the CO2 fractional abundance to ~ 3 x 10~ 5 is discussed below. 

To simplify the calculation we assume that the average temperature is below 30 K and 
that the hydrogen density is constant throughout the envelope. At this low temperature 
thermal desorption of CO2 is negligible and therefore the equilibrium gas phase abundance 
of CO2 is dependent only on the UV photodesorption and the freeze-out rates. We also 
make the approximation that the total CO2 abundance is constant throughout the envelope, 
since observations show that CO2 forms in ices at the edges of clouds and does not increase 
deeper into the cloud (jPontoppidanl 120061 ) . The UV field is composed of the interstellar 



radiation field of 10 8 photons cm~ 2 s , which is attenuated with Ay, and the UV photons 
produced inside the cloud by cosmic rays. For a typical galactic cosmic ray flux, the resulting 
UV photon flux is of the order of 10 4 photons cm -2 s _1 with a factor of 3 uncertainty 



(Shenetal 



20041 ). Using our derived photodesorption yields and the estimated freeze-out 



rate of CO2 we calculate the steady-state gas abundance of CO2 as a function of Ay for 
small and large grains. 
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The photodesorption rates of CO2 molecules from grain surfaces in molecules s _1 due 
to external and cosmic ray induced UV photons, respectively, is described by 

-RpD-ISRF = -^ISRF-FUV e_7Av ^pd ( na gr) ! (2) 
-RpD-CR = ^CR-FUV^pd ( 7ra gr) ; (3) 

F pd = 0.0012(1 - e- x/2 - 9 ), (4) 

x = n c02 -ico/ (1 x 10 15 x na 2 gr n gT ) , (5) 
and 

"co 2 -icc = 5x 10~ 6 n H - n C o 2 -gas, (6) 

where /pd-isrf is the strength of the external irradiation field with energies 6-13.6 eV 
and ipD-FUV is the strength of the UV field due to cosmic rays. 7 is a measure of UV 
extinction relative to v isual extinction, which is ^2 for small interst ellar grains and <0.6 



for grains of a few fjwi (jRoberge et al.lll991tlvan Dishoeck et al.ll2006l ). and a gr is the grain 



radius. V p d is the experimentally determined CO2 photodesorption yield for temperatures 
below 30 K (Eq. 2]). x is the ice thickness in monolayers, which is defined in terms of the 
CO2 ice abundance (n.co 2 -icc in cm -3 ), grain surface (7ra 2 r n gr in cm 2 cm -3 ) and amount 
of molecules per monolayer (1 x 10 15 cm -2 ) in Eq. [5l Equation [5] states the relationship 
between the gas and ice phase CO2 abundances (nco 2 -gas and «co 2 -icc) when the total 
CO2 abundance is 5 x 10~ 6 nn and nn is set to 1 x 10 4 cm~ 3 . 

The accretion rate of CO2 molecules i? acc is a product of the molecular velocity, the 
grain surface, the sticking coefficients and the gas phase abundance of CO2 according to 



R acc = -4.57 x 10 4 ) ' (vra 2 ) Sn C o 2 -gas. (7) 

\mc02J 

In Eq.[7]the gas temperature T is set to 10 K, mco 2 is the CO2 mass in atomic m ass units, 



2006). 



and S is the sticking coefficient, which is assumed to be 1 (Bis schop et al. 

At steady-state the total photodesorption rate (i.e. the sum of Eqs. [5]and[3]) is equal to 
the accretion rate. Figure rjj] shows the resulting steady-state CO2 gas phase abundance 
as a function of Ay for small classical 0.1 /xm grains and after grain growth to a few /xm. 
At Ay <3 mag, CO2 gas is photodissociated and the model is not valid there. Deep into 
the envelope, the gaseous CO2 abundance due to cosmic ray induced photodesorption is 
~ 2 x 10~ 8 , which is about an order of magnitude lower than the observed CO2 abundance. 
For small grains the external UV light also probably does not penetrate deep enough into 
the envelope to increase the average abundance significantly. On average about 1% of the 
total CO2 abundance is kept in the gas phase through photodesorption. If the CO2 ice 
abundance is increased to ~ 3 x 10~ 5 , the CO2 gas abundance increases by an order of 
magnitude at low extinctions. Deep into the cloud the abundance only increases by a factor 
of two, since the photodesorption rate only depends on the ice thickness up to ~10 ML. 
In cloud and envelope material dominated by small grains, photodesorption probably does 
not explain gaseous CO2 abundances of 5-10%, as observed toward L 1527 IRS, unless 




Fig. 14. The CO2 gas phase abundance with respect to the total hydrogen (H+H2) column 
density assuming a total (gas 4- ice) CO2 abundance of 5 x 10~ 6 . The two tracks assume 
classical 0.1 /jm grains and grain growth up to a few /nm, respectively. 



there is a str ong internal UV so urce and the UV is scattered into the cavities created by 
the outflow (ppaans et al.l Il995f) where it can photodesorb material in the cavity walls. 
Indeed, LI 527 IRS is well known for its pr ominent outflow and X shaped cavity wall o n 



scales comparable to the IRAM 30m beam ([MacLeod et al 



1994; 



Hogerheiide et al 



1998). 



In disks, or in general where grain growth has occurred, the picture changes dra- 
matically. Millimeter observations of oute r disks show grain growth up to mm size 



(Rodmann et al 



2006 



Lommen et al 



20071 ). and the UV irradiation may then penetrate 
deep into the disk releasing a high fraction of the CO2 ice into the gas phase. At an Ay 
of 10 mag the CO2 fractional abundance is now ~ 1 x 10~ 6 or 20% of the total CO2 
abundance. Similarly to the case with small grains, increasing the CO2 ice abundance only 
increases the CO2 gas abundance substantially at low extinctions. 

This simple model thus shows that the experimentally determined CO2 photodesorption 
yield is high enough to release large amounts of CO2 ice into the gas phase if moderate 
grain growth has occurred. Even with small grains and no internal UV source except that 
produced by cosmic rays, photodesorption may keep up to ~1% of the total CO2 abundance 
in the gas phase deep into cloud cores. This model stresses, once again, that photodesorption 
of ices needs to be taken into account when modeling gas-grain interactions. 



5. Conclusions 

1. The CO photodesorption yield is temperature dependent between 15 and 27 K, which 
is described empirically by 2.7 x 10~ 3 — (T — 15) x 1.7 x 10~ 4 molecules photon -1 . The 
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anti-correlation between yield and temperature is probably due to ice re-structuring into 
a more compact configuration - the observed linearity may be coincidental, however. 
For most astrophysical applications the yield measured at 15 K is appropriate to use. 

2. The CO photodesorption is initially reduced by more than 80% when the CO ice is 
covered by 1 ML of N2 ice and decreases with UV fluence when mixed with N2 due to 
surface build-up of N2 ice, confirming that CO only desorbs from the ice surface. 

3. N2 co-desorbs with CO at 16 K in an ice mixture and in a layered ice with a yield of 
3 x 1CP 4 molecules photon -1 . 

4. A CO2 photodesorption event starts with the photodissociation of a CO2 molecule into 
CO and O. The fragments either desorb directly or react and recombine to form CO2 
and CO3 before desorbing. The CO3 yield is however less than 1% and the two main 
desorption products are CO and C02- 

5. The CO2 photodesorption yield is thickness dependent at all temperatures between 18 
and 60 K. At 18-30 K the yield is well described by 1.2 x 10~ 3 x (1 - e -*/ 2 - 9 ) + 1.1 x 
10~ 3 x (1 - e-*/ 4 - 6 ), and at 40-60 K by 2.2 x 10" 3 x (1 - e - x / 5 - 8 ) + 0.22 x 10~ 3 x x 
molecules photon -1 , where x is the ice thickness in monolayers. The first part in each 
yield equation is due to desorbing CO2 molecules and the second part to desorbing CO 
molecules. 

6. The thickness dependence of CO2 photodesorption is understood from a mean-frce-path 
perspective, where the different excited fragments travel a different average distance 
through the ice before being stopped. At higher temperatures, this mean free path 
increases due to increased mobility of molecules in the ice. 

7. A simple model of an envelope using the observed CO2 abundance in L 1527 IRS shows 
that CO2 photodesorption can maintain CO2 fractional abundances up to 1 x 10~ 6 in 
the gas phase at Ay ~ 10 mag after moderate grain growth and 2 — 3 x 10~ 8 using 
small ISM grains. At lower extinctions the photodesorption is higher due to the external 
irradiation field and a high fraction of the total CO2 ice abundance is maintained in 
the gas phase. 
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